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Abstract: We present a simple derivation of the stochastic equation obeyed by the 
density function for a system of Langevin processes interacting via a pairwise potential. 
The resulting equation is considerably different from the phenomenological equations 
usually used to describe the dynamics of non conserved (Model A) and conserved 
(Model B) particle systems. The major feature is that the spatial white noise for this 
system appears not additively but multiplicatively. This simply expresses the fact that 
the density cannot fluctuate in regions devoid of particles. The steady state for the 
density function may however still be recovered formally as a functional integral over 
the coursed grained free energy of the system as in Models A and B. 
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Recently there has been a tentative attempt [T|] to analyse the liquid-glass transition 
within the context of the replica approach to disordered systems. A program of imple- 
menting replica techniques from disordered systems has had some success in capturing 
the properties of systems without disorder but which have many metastable states 
|2|, |3], |]]. As the glass transition is essentially dynamical it is useful to have a simple 
dynamical model for these types of problems. In the literature on liquids there has 
been a considerable effort to analyse the dynamics of liquid systems. In this letter 
we present a simple microscopic derivation for the evolution of the of the density of a 
system of Langevin processes interacting via pairwise potentials. Rather by working 
with Fokker-Planck equations in large coordinate spaces and then using projection op- 
erator techniques to obtain a hierarchy of equations for the n point density correlation 
functions we simply derive a closed functional Langevin equation for the evolution of 
the density. Using this it is the trivial to derive the dynamical BBGKY (e.g. see || [5| 
and references therein) hierarchy of equations. Interestingly the naive course grained 
free energy of the system appears quite naturally. However the equation is not amongst 
the class of equations used often to describe the approach to equilibrium of systems 
specified by their free energy as a function of some scalar order parameter |?|. |8| . 
We commence by considering the dynamics of a system of particles in a thermal white 
noise heat bath. The particles interact via a pairwise potential V(x), this could be a 
Coulomb type of interaction, a Leonard-Jones type interaction or simply a hard sphere 
repulsion. Each particle experiences a thermal white noise and moves under the force 
generated from the potential due to its neighbours. Consequently the ith particle obeys 
the Langevin equation 

HX (t) N 

= m(t) - 5yv(x t (t) - x 3 (t)), (i) 

to ease notation we have assumed a potential such that V^(0) = 0. The noise is 
uncorrelated in time and the noise acting on a single particle is not correlated with 
the noise acting on the others. The components of the noise are also taken to be 
uncorrelated, hence: 

(r ] mV-(t')) = 2T5 l ^5(t-t>). (2) 
Our strategy will be to consider the evolution of the density function for a single particle 

Pi(x,t)=6(Xi(t)-x), (3) 

we shall then demonstrate how one may write a closed Langevin equation for the global 
density 

N 

POM) = 5> i (x,t). (4) 

The derivation follows a well known argument. Consider an arbitrary function / defined 
on the coordinate space of the system. Using the definition of the density it is a 
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tautology that 

f(X i (t))= I dx Pl (x,t)f(x). (5) 



Expanding the stochastic differential equation using the Ito calculus (for example see 
H) over the next time step St one obtains 



dt 



J dx Pi {x,t) (vf(x)-TH{t)-Vf{x) ■ (fyv(x - Xjit))) + TV 2 f{x) 



(6) 

Rearranging the above integral by integration by parts we obtain 



df(Xj 
dt 



J dxf(x) ■ (piix,^^)) + V ■ Pi (x,y) \Y,VV(x - X^tyj + TV 2 Pi (x,t)} . 

(7) 



However from (y) we may also deduce 



df (x i )_ r dx d_Pd^ f{x) 



dt J dt 

Comparing equations ([7|) and (§) we find (using the fact that / is an arbitrary function) 
that 



dpi(x,t) 



N 

()f V-(p,(x,t)^(t)) + V- [(h(x,t)Y l VV(x-X j (t))) +T\7 2 Pt (x,t). (9) 



We emphasise that this argument is standard and the only subtlety is that we have 
not carried out any thermal averaging at this point. Summing equation (|9]) over the % 
and using the definition of the density p we obtain 



dp(x,t) 
dt 



" E V • (Pi& tMt)) + V ■ (p(x, t) f dyp(y, t)VV{x - y)) + TV 2 p(x, t). 

(10) 

This equation is almost a closed equation for p, the problem is that the noise term 
appears to contain too much information about the individual pi and hence, in the 
current form, it is not a Markovian equation for the evolution of the global density. 
We define this noise term to be 

N 

t(x,t) = -Y,V-{Vi(t)Pi(x,t)), (11) 

i=i 

it is clearly still Gaussian and the correlation function is given by 

N 

(£{x, 0) = 2TS(t - t>) ]T V x ■ V y ( Pi (x, t) Pi (y, t)) . (12) 

i=l 
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However the above term simplifies using a trivial property of the Dirac delta function, 
explicitly 

Pi(x,t)pi(y,t) = 5(x - y)pi(x,t) = S(x - y)pi(y,t). (13) 
Using this one finds 

(C(x, t)£(y, t')) = 2T5(t - t')V x ■ V y (S(x - y)p(x, t)) . (14) 

At this point we shall consider redefining the white noise by introducing a global noise 
field 

^'(x,t)=V-( V (x,t)pkx,t)), (15) 
where i] is a global uncorrelated white noise field such that 

(rf(x, t)T} v {y, 0) = 2T5(t - t')5^5(x - y). (16) 

It is easy to see that both Gaussian noises £ and £' have the same correlation functions 
and are therefore statistically identical. We may therefore write equation ( |i~0|) as 

?P^ = V-( V (x,t)p^x,t))+V-(p(x,t) J dyp(y,t)VV(x-y^+TV 2 p(x,t). (17) 



Equation ( |TTD is clearly a closed Langevin equation for the evolution of the particle 
density of the system. The crucial point is that on insisting on knowing the global 
property of the white noise field we have gained sufficient information to write down 
a Markovian evolution equation for p. If one were to consider the course grained free 
energy for the system one would be lead to write 

F = — f dxdy p(x)V \x — y) p{y) + T / dxp(x) log p(x) , (18) 



up to constant terms in the case of a fixed particle number. Using this as a definition 
of F it is easy to see that ( |TTD may be written as 

= V - Uar,t)V^y| p(a!lt )J +V - (p(x,t)p*(x,tj) . (19) 

Two things are immediately obvious from equation (IDp, firstly particle number is 
conserved and secondly in regions where the system is void of particles the density 
cannot change and hence stays zero. The above equation is in fact a special case 
of an equation derived semi-phenomenologically by Kawasaki some 20 years ago in 
the context of non-polymeric incompressible binary fluids [10|. As is usual in such 
derivations, the free energy is taken to be of the Ginsburg-Landau type directly, if 
one sets the velocity field to zero one obtains formally the same equation as (|17|). The 
derivation in the case of zero velocity field is brief enough to be included here for the 
sake of completeness. One assume that the density obeys an equation of the form 

§j = -V-j + £0M), (20) 
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where j is the current and £ is the divergence of a random force uncorrelated in time. 
In the course grained microscopic approach j is simply given by 

i = -pvf , (21) 

the fluctuation dissipation theorem (see for example the discussion in [|TT[]) then requires 



that the noise £ must have exactly the same correlation function as the noise in flT7|) 
hence the result. 

Here we briefly recall the two standard models used in dynamical analysis of phase 
ordering kinetics and other areas [|7|, |§. The basic approach here has to been to write 
down a course grained free energy F for the system and to write for Model A (non 
conserved dynamics) 

<*> 

and for Model B (conserved dynamics) 

|.V^ + V.,(M). (23) 

Here p is not necessarily the density and is usually some scalar order parameter for the 
problem, for example in a binary alloy it would correspond to the difference between 
the densities of the two particle species and in a spin model it would be the local 
magnetization. In Model B the noise is vector white noise of the type introduced above 
and in Model A it is the scalar counterpart. The models can lead to an equilibrium 
distribution for p of the type that has as its free energy the F one started with, i.e. 
(formally) the equilibrium functional for p is given as 

P(p) = ^ex P (-F[p]/T) (24) 

with Z a normalizing factor. One may verify that the functional Fokker Planck equation 
for equation ( |17D also has this as a steady state solution. In this case it is interesting 
that the course grained free energy has appeared directly and in addition without any 
appeal to there being a large number of particles in the system, we recall that the term 
/ dxp(x) log p(x) normally occurs after a course graining and then a use of Stirling's 
formula. In this formulation we have avoided this route, presumably it is the direct use 
of the Langevin process at the level of microscopic dynamics that has smoothed out 
the statistics. However the functional representation of P[p] does require additional 
reference to the underlying physics when it comes to performing what is an a priori 
badly defined functional integral. In the case of both Model A and Model B (and 
indeed for the statics) the usual tactic is to resort to an expansion about some fixed 
background density up to fourth order in the fluctuations, this then yields a dynamical 
version of Ginsburg-Landau theory. However we note that, even at the level of this 
approximation, equation ( ]T7| ) does not become equivalent to Models A or B (at this 
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level of approximation one would still be forced to keep terms of third order in the 
terms multiplying the white noise). 

Once we have the full Langevin equation for the density evolution, we may use the 
standard techniques of stochastic calculus (applied now to functions) in order to obtain 
equations for the evolution of functions of the density. Here we shall consider the time 
dependent correlation function at equal times. We define 

C(x,y,t) = (p(x,t)p(y,t)). (25) 

It is straightforward to see that C(x,y,t) obeys 

dC{x,y,t) 
dt 

- TV 2 6(x - y) (p(x, t)) +2(V ■ t)p(x, t) J dzp(z, t)VV(x-z)^} + 2TV 2 C(x, y, t) , 

(26) 

where we have use the symmetry between x and y. Unfortunately the equation (^6j) is 
not a closed equation for C. In general the above may be used to generate a whole hier- 
archy of equations for the n point functions of the field p, in fact this is the dynamical 
version of the BBGKY hierarchy used in the static problem. 

In conclusion we have presented a very simple derivation of the Langevin equation of the 
density describing the dynamics of a gas of interacting Langevin processes. In addition 
this formalism permits the construction of the BBGKY hierarchy for the dynamical 
correlation functions in straightforward manner. Rather than being a phenomenological 
model for the dynamics based around a known course grained free energy, it is directly 
defined from an albeit simple microscopic model. It has the merit of corresponding 
to a (at least in the sense of numerical simulations) realisable microscopic dynamics. 
One rather appealing feature is that the naive course grained free energy of the system 
appears naturally as the effective potential for the problem. 
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G Parisi, and support from EU grant CHRX-CT93-0411. 
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